Sensitivity of Wave Field Evolution and Manifold Stability in Chaotic Systems 



Nicholas R. Cerruti and Steven Tomsovic 

Department of Physics, Washington State University, Pullman, WA 99164-2814 

(Dated: February 8, 2008) 

The sensitivity of a wave field's evolution to small perturbations is of fundamental interest. For 
chaotic systems, there are two distinct regimes of either exponential or Gaussian overlap decay 
in time. We develop a semiclassical approach for understanding both regimes, and give a simple 
expression for the crossover time between the regimes. The wave field's evolution is considerably 
more stable than the exponential instability of chaotic trajectories seem to suggest. The resolution 
" l"i of this paradox lies in the collective behavior of the appropriate set of trajectories. Results are given 

' for the standard map. 
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^ I There is a multitude of physical systems in which the quantum evolution of a particle or that of a wave field 
. interacts with a time-varying, complicated, essentially-unknowable medium. Well known examples include acoustic 
waves propagating long ranges through the ocean's sound channel Q, starlight passing through our atmosphere 
and electrons diffusing through a metal [^. In a similar vein, wave/quantum evolution problems whose underlying 
ray/classical analogs are simple chaotic systems, may possess parameters which are tunable or time varying as well. 
l-^ . In all of these systems, understanding the effects that system changes have on the dynamics is of great interest. 
' In this letter, we investigate the problem of the overlap decay of two initially identical wave fields propagating in 
slightly different simple, chaotic systems and using semiclassical dynamics derive the functional dependence of the 
decay process on time, Planck's constant, and perturbation strength. 
P5 . Recently, the overlap between time-evolved perturbed and unperturbed states has been studied in many contexts 
and guises. The overlap has been suggested as a measure for the stability of quantum motion It was identified 
I . with the Loschmidt echo and studied in the context of polarization echoes in nuclear magnetic resonance ^ where 
^ ' a cross-over between Gaussian and exponential behavior of the overlap was observed. In quantum computing, it is 
• equivalent to the fidelity which measures the loss of phase coherence |7|. The fidelity was shown to have different 
functional decay forms in the chaotic and integrable limits |^ . 

Over the past ten years, it has been shown that semiclassical theories of wave mechanics based on rays/trajectories 
' are capable of quantitatively reproducing evolving wave fields in chaotic dynamical systems for significant propagation 
y—i , times during which the chaos can become extremely well developed |lO| . Extensions of these semiclassical methods 
should therefore provide an accurate theoretical approach for understanding overlap decay behavior to at least the 
, same significant times. We therefore couple a semiclassical approach with classical and quantum perturbation theory 

• '-H ' as appropriate for the small system changes we are considering, and methods from random matrix theory, which are 
^ ■ known to apply to fully chaotic systems. 

^ ' A number of paradoxes arise in considering the semiclassical limit. For example, in a classically chaotic system a 

• ^ , slight perturbation of the initial conditions causes an exponential divergence of the rays and their action differences; 
' the exponent characterizing this divergence is the Lyapunov exponent. In addition, perturbing the Hamiltonian, but 

using the same initial conditions, also causes the rays to exponentially diverge with the same Lyapunov exponent as 
for a change in the initial conditions. Wave fields behave very differently under perturbations than rays. The same 
wave field propagated in two slightly different systems has an overlap which decays in time on a scale proportional to 
the square of the strength of the system changes. On the other hand, two wave fields starting from different initial 
conditions in the same system have an overlap which is constant in time. Those paradoxes that exist must have 
resolutions or the semiclassical theories could not be such faithful approximations. Their resolutions hold the promise 
of revealing new features of either classical or wave mechanics. We discuss one example at the end of this letter. 

Consider an evolving quantum wave function \a(t)x) where A represents parameters that define the system. The 
overlap decay is defined as: 

C{e;t)^Cx{e;t) = \{a\Ul_^^,{t)Ux+,/2{t)\a)\' (1) 

where e serves as the small parameter; the use of ±e/2 reduces next-to-leading perturbative corrections. C{e;t) can 
be conceptualized equally well the squared overlap of an initial state \a) propagated forward in time with two slightly 
different Hamiltonians or as |a) propagated forward in time with one Hamiltonian and then backward in time with a 



o 

00 



2 



slightly different one. In the limit that the system is strongly chaotic, the statistical properties of C\{e;t) would be 
stationary rendering A largely a superfluous subscript. 

Assuming the Hamiltonian is differentiable with respect to A, let V = dH{\) / dX. The simplest approximate 
expression for C(e; t) follows by keeping only the leading quantum perturbation corrections to the energies since they 
are divided by Ti; this is clearly insufficient to describe fully the overlap decay behavior, but it does describe the first 
limiting regime. The resulting expression is 
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where ban = {a{Q)\nx) and Km = ("-a|V^|"-a)- 

In the strong chaotic limit, Eq. (|^), can be evaluated using statistical arguments. Random matrix theory asserts 
that the eigenvalues and eigenfunctions are uncorrelated, and thus the amplitudes and phases statistically decouple; 
semiclassical arguments have also been developed that support this approach. The phases are Gaussian random 
distributed and are proportional to the derivative of the eigenenergies with respect to the parameter, commonly 
called the level velocities. Thus, the sum can be replaced by an integral over the Gaussian density, which gives 

C(e;0 «exp(-eV|tV?i') (3) 

where <7^ is the variance of the level velocities. In the small e and strongly chaotic limit, the overlap decay has a 
Gaussian form with respect to time. 

The level velocity variance is not a free parameter, and has been related to V through a semiclassical approach |ll| , 
111; it is 
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where g is a symmetry index, d is the mean density of states, and K{E) is the classical action diffusion constant given 
by [pl 



K{E) = / {V{p{0), q(0); X)V{p{t), q(t); A)), 
Jo 



dt 



(5) 



The averaging is defined over the primitive periodic orbits of very large period. 

Equations (|^^ turn out to describe a very restricted range of e centered at zero. A semiclassical approach without 
the quantum perturbation limitations can be developed that is valid over a much broader range. The semiclassical 
construction of an evolving wave function begins with the propagator 
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The phase is specified by the time integral of the Lagrangian Wj{q, q'; t) and an index based on the properties of the 
conjugate points (like focal points), Vj. 

The overlap decay in terms of the propagator is 



C(e;t) 
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For simplicity, we have chosen the initial states to be point sources, i.e. {q\a) = 5{q — go)i but the following theory is 
generalizable to any initial state. 

The most important contribution to the overlap decay is each orbit's action change due to its division by %. Thus, 
changes in the amplitudes are ignored, so 
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FIG. 1: Examples of Gaussian, exponential and cross-over decay behaviors. The quantum kicked rotor curves are solid, and 
the theoretical curves dashed. The parameters are N = 350, A = 18, g = 4, K{E) = (1 + 2J2(A))/4(27r)'', and in (a) e = lO"*, 
in (b) 5 * 10""^, and in (c) 1.5 * 10~^. Note that the equivalent expression of Eq. for maps is a|; = 2gK{E)/N. 



The AWn{q, qo',t) are the action differences of two rays that begin at go and end at position (7 in a time t and are 
continuously deformable into each other; bifurcations are neglected. First order classical perturbation theory gives 

Aw^4.-,.;t) = e{ '^"^^;^^^'^ c^.- (9) 

where the amplitude in Eq. and the Lagrangian are evaluated along the orbits of H{\). In general, stationary 
phase integration cannot be performed on Eq. (^) because the action differences are less than fi. 

In strongly chaotic systems or random media, the classical action differences for long orbits obey a central limit 
theorem. The contributions of the derivative of the Lagrangian along parts of an orbit separated widely in time are 
uncorrelated. Furthermore, each orbit has a unique history in which it visits the available phase space. Thus, the 
changes in action will be Gaussian distributed for sufficiently long times. The integral over q can be replaced by a 
Gaussian probability integral over the action change which gives 

C(e;i) «exp(~a^/?i2) (10) 

CT^ is the variance of the action differences and is related to the same classical diffusion constant K{E) as are the 
level velocities 



= 2e^K{E)t 



(11) 
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FIG. 2: The same two initial manifolds {qo — 0.6) propagated four time steps through two slightly different media in the 
quantum kicked rotor where Ai = 6.0 and A2 = 6.02. The left panel is on the torus and the right panel is unfolded in phase 
space. 



In contrast to the quantum perturbative argument giving a quadratic i-dependence in the exponential, here the 
argument is linear. 

The exponential or Gaussian dependence dominates over the time range in which the argument of the exponential 
takes on the lesser value. The cross-over between the two regimes occurs when their arguments are equal. Simple 
algebra gives the cross-over time scale t* — hd/g, which is essentially the Heisenberg time, th- Exponential decay 
dominates if the decay is completed by th , and Gaussian decay dominates if little decay has occurred by th ; otherwise 
the decay has components of both behaviors. In terms of the parameters, the question is 



e2 > ^ = -J^ (12) 

< 2K{E)th indKiE) ^ ' 



This cross-over fits naturally with the approximations made in deriving the exponential and Gaussian behaviors. 
The classical action correlations necessary for a semiclassical theory to "know" the spectrum is discrete were not 
included whereas the simplified quantum perturbation theory left out coupling with neighboring levels which would 
be completely dephased by th- 

We use the standard map, a paradigm of chaos for large kicking strength, to illustrate our main results. The 
classical map is defined by 

Pi+i = _Pi - (A/27r) sin(27rqi) mod(l) 

qi+i = qi + Pi+i mod(l) (13) 
For A greater than 6 or so, the system is strongly chaotic. The quantized propagator with N discrete levels is given 

by 

(n'|[/|n) = -^exp[iTr{n-nf/N] 

X exp ( j— cos[27r(n' + a)/Ar] j (14) 

where n,n' — 0, . . . , N — 1 and a is a phase term which we set equal to zero. The effective Planck constant is given 
by /i = l/N. 

In Fig. ^ we show three example overlap decays. The upper panel illustrates an example of Gaussian decay, the 
middle panel illustrates exponential decay, and the bottom panel illustrates a case where the cross-over between the 
two cases can be seen. The theoretical curves drawn for comparison have the analytic value of K{E) |ll] inserted, 
and so involve no parameter fitting whatsoever. The agreement is excellent. 

At first consideration, it may seem hopeless to apply a first order classical perturbation theory to a chaotic system. 
The action difference cannot grow faster than the time, whereas the classical action of a trajectory associated with 
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FIG. 3: The action changes along the manifold. The solid line in the the left panel is the exact action differences for the two 
manifolds in Fig. 2| and the dashed line is the perturbation result. The right panel is a histogram of the action changes for 
eight time steps of the initial manifolds in the right panel and the solid line is a Gaussian fit to the histogram. 



a particular initial condition deviates exponentially. At best, the first order approximation should be valid no longer 
than a logarithmically short time scale, r* = — i Ine, and the theory given above should not work. The resolution rests 
in considering the difference between the behavior of a single trajectory and the collection of trajectories underlying 
the construction of a wave field. In time-dependent WKB theory, the collection is a Lagrangian manifold. In our 
example here with a position state initial condition, the initial manifold is just the set of initial conditions containing 
all momenta and q = qq. To propagate the state, the manifold is propagated. Likewise to distinguish the propagation 
of the same initial wave function for slightly different initial conditions entails propagating the manifold under the 
two slightly different dynamics. This is illustrated in the left panel of Fig. |^ for an initial g-manifold. The two 
manifolds are so close as to be nearly indistinguishable. In the bottom panel where the exact same propagated 
manifolds are depicted except with the phase space torus unfolded into a plane, the -l-'s mark the final point of a 
single initial condition. The final point deviates an entire unit cell away, and yet has only essentially slid along the 
original manifold. The key feature of the slowly varying manifold is related to the changing locations of the folding 
points, which generally are associated with caustics. The manifold has a strong structural stability that still allows 
for exponential sensitivity to perturbations for its individual members. 

Thus, the construction of an evolving wave function reflects the stability of the manifold, not the individual trajec- 
tories. For example, in the Green function, the condition of starting at q' and ending at q naturally selects trajectories, 
continuously deformable into each other, whose action differences under perturbation are minimal. The first order 
classical perturbation theory accounts for this effect. In Fig. |[ the action difference is plotted in the left panel as a 
function of position along the manifold. The first order classical perturbation prediction is superposed, and it is im- 
possible to distinguish them; it is extremely accurate. It turns out that the extrema of the action changes correspond 
to the points where the two manifolds intersect. In the right panel, we see that the fluctuations are such as to give a 
good approximation to a Gaussian density. In the limit of longer orbits, the results would get better and better. 

In conclusion, the overlap of a wave function evolving under two slightly different chaotic dynamics decays either 
exponentially or in a Gaussian fashion. The decay scales and cross-over behavior can be derived semiclassically. We 
have not shown it here, but deviations due to not being 'chaotic enough' or other statistical shortcomings are captured 
by the semiclassical theory. However, the correlations and deviations are buried in the orbit sums, and not easily 
extracted other than by numerical summation. An interesting consequence of classical mechanics is that manifolds 
are much more stable against perturbation than the trajectories of which they are comprised. First order classical 
perturbation theory yields an excellent approximation to the changes of the actions, and a statistical theory of these 
action variations leads to a very good agreement of the overlaps. The stability of wave fields is connected to the 
manifolds and the associated phases, and not the individual trajectories. 

We gratefully acknowledge important discussions with M. A. Wolfson and T. Prosen, and support from the Office of 
Naval Research under Grant No. NSF-PHY-0098027 and the National Science Foundation under Grant No. N00014- 
98-1-0079. 

Note: After this work was completed, we discovered Ref. [|4| which also examines the overlap decay. 
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